Lifetimes of electrons in the Shockley surface state band of Ag(lll) 
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We present a theoretical many-body analysis of the electron-electron (e-e) inelastic damping rate 
F of electron-like excitations in the Shockley surface state band of Ag(lll). It takes into account 
ab-initio band structures for both bulk and surface states. V is found to increase more rapidly as a 
function of surface state energy E than previously reported, thus leading to an improved agreement 
with experimental data. 
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I. INTRODUCTION 

The quantitative understanding of the dynamics of 
electronic excitations in the Shockley state band at the 
(111) face of Ag has been the subject of a considerable 
number of research reports (for a review see Ref. pj). 
Recent experimental studies used photoelectron spec- 
troscopy and the scanning tunnelling microscope (STM) 

Theoretically, many-bod y li fetime calculations using 
the GW approximation [3, [H| have clarified the role of 
electron-electron (e-e) interactions as well as the relative 
importance of intraband and interband scattering and 
screening. The damping rate due to electron-phonon 
scattering is also quantitatively understood How- 
ever, comparing theoretical with STM-derived experi- 
mental damping rates of electronic excitations with en- 
ergies above the Fermi level Ep show that the experi- 
mental decay rates increase more rapidly with increasing 
energy E. Recently a similar discrepancy between im- 
age state lifetimes measured by means of the STM and 
two-photon-photoemission spectroscopy was lifted [l2| . 
It was shown that the applied electric field during the 
STM- measurement causes an increase in both the effi- 
ciency of the image state decay channels as well as their 
number. In the light of this work, Becker et al. theoreti- 
cally investigated the impact of the STM-induced electric 
field on the inelastic e-e damping rate of the electron- 
and hole-like excitations in the Shockley surface state 
band at Ag(lll) |ljj. Their results indicate that under 
typical tunnelling conditions the STM does not signifi- 
cantly alter the surface state wave function and so pre- 
vious STM-derived Shockley state lifetimes need not be 
corrected. Nevertheless, while experimental and calcu- 
lated lifetimes agree quite well over a range of energies, a 
quantitative difference remains for unoccupied states of 
energies E > Ep + 300 meV. The origin of this differ- 
ent energy dependence, which is still an open question, 
is addressed in this report. 



II. CALCULATION DETAILS 

The calculations presented here are based upon the ap- 
proach developed by Chulkov et al. |14| , and used widely 
in calculations of surface state dynamics [l| . Within this 
approach the damping rate or inverse lifetime of an exci- 
tation in the state ip(r) with energy E is obtained from 
the expectation value of the imaginary part of the elec- 
tron self-energy, E(r, r';E): 



r = r- 1 = 



= -2 Jdr Jdr' ijj*(r)ImT,(r,r'; E)tp(r'). 



(1) 



Unless stated explicitly, atomic units are used throughout 
the text, i.e., e 2 = h = m e = 1. 

In the GW approximation ^| of many-body theory the 
imaginary part of the self energy is calculated to first 
order in terms of the screened Coulomb interaction W 
and the Green function G 

Im£(r,r';e) = -- f de' ImG(r, r'; e')lmW(r, r'; e-e'). 

I" JE F 

. (2) 

Here we apply the widely used zeroth order approxima- 
tion to the Green function in the spectral representation 



^(r)^(r') 
e - E t + iS ' 



(3) 
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where the tpi(r) are one-electron eigenfunctions with en- 
ergies Ei and 5 a positive infinitesimal. The screened 
Coulomb interaction is evaluated in the random phase 
approximation (RPA) 

W(r, r';w) = v(r — r') + J dri J dr 2 v(r — i^) 

x X °(ri,r 2 ;u;)Ty(r 2 ,r'; W ) (4) 

where v is the bare Coulomb interaction and x° is the 
density-density response function of the non-interacting 
electron system. An explicit expression for x° can be 
found in Ref. [Ttij . 
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FIG. 1: (Color online) Ag(lll) surface band structure. 
Shaded area and solid line are the projected bulk and surface 
state dispersion from ab-initio calculations which we use for 
calculating the e-e decay rate. The dashed lines are parabolic 
dispersions with effective masses m* = 0.4 (surface state 
band) and 0.25 (band edge). 



The single particle states have been obtained by solv- 
ing the Schrodinger equation using a one-dimensional 
pseudopotential varying in the direction perpendicular 
to the surface 01- The pseudopotential does not de- 
scribe the d-electrons of the substrate, which are any- 
ways too low in energy to play a significant role as final 
states for the decay, but their contribution to the screen- 
ing is included via a polarisable background . The 
electron-phonon contribution to the linewidth is taken to 
be constant r o _ p h = 3.7 meV for energies in excess of 
20 meV above the Fermi level [TT| . 

Previous calculations of the decay rate of the Ag(lll) 
Shockley state have assumed parabolic dispersion with 
effective masses of m* = 0.4 and m* = 0.25 for the in- 
trinsic surface state and the lower band edge respectively 
However, over the extended energy range of interest 
here, this approximation needs to be improved as evident 
from Fig. 2] Figure ^ shows the parabolic dispersions 
(dashed lines) for the surface state band and the band 
edge with effective masses as mentioned above. The grey 
shaded area and the solid line are the projected bulk and 
surface state dispersion, respectively, that are found from 
ab-initio calculations. The parabolic dispersion underes- 
timates the gap between surface state band and band 
edge. We have therefore used the ab-initio dispersions 
for our lifetime calculations. As previously noted |sL ITflj 
there are important changes in the shape of the surface 
state wave function with k, the surface wave vector, and 
we take these into account by recalculating the wave func- 
tion for different k with the pseudopotential parameters 
changed to take into account the appropriate ab-initio 
band edges and surface state energy. 



FIG. 2: (Color online) Comparison of the total Shockley sur- 
face state decay rates r+r e _ p h for electrons at Ag(lll) calcu- 
lated using ab-initio dispersion relations (solid line) and previ- 
ous published decay rates using parabolic dispersion (dashed 
line) from Ref. fil. The black squares represent experimental 
data from Ref. [g, black triangles from Ref. and black cir- 
cles from Ref. The experimental data from Refs. 0,0 have 
been multiplied by a factor 2 to compensate for an erroneous 
analysis in that work 20]. 

III. RESULTS AND DISCUSSION 

We focus on the energy range between and 1 eV above 
the Fermi level where the surface state is energetically 
separated from the bulk states and thus is well defined. 
This energy range is of interest as according to Ref. Q 
there is a cross-over in the primary decay channel from 
intraband to interband dominated decay. Our calculated 
results for the total damping rate r + r c _ p h are shown in 
Fig. [2] as a solid line. For energies below E — Ep ss 0.1 eV 
we find that our results for the inelastic e-e damping 
rate of excitations in the Shockley state band agree well 
with the inverse lifetime estimation from Echenique et 
al. (dashed line) [l| using similar methods but with an 
effective mass approximation. This agreement is due to 
the fact that the effective mass approximation gives a 
very good description of the projected ab-inito surface 
band structure for energies with E — Ep < 0.1 eV as 
evident from Fig. ^ 

For larger energies we find that the e-e damping rate 
as a function of energy increases more rapidly with the 
ab-initio dispersion than calculated with parabolic dis- 
persions. This result can be understood as follows: For 
energies E — Ep > 0.1 eV the parabolic dispersions un- 
derestimate the energy gap between surface state band 
and the bulk band edge. Consequently, the ab-initio dis- 
persions reveal a reduced penetration of the surface state 
wave function into the interior of the crystal and an in- 
creased probability amplitude in the proximity of the 
metal surface. This modification leads to two compet- 
ing effects. Since screening is reduced at the surface [E| 
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intraband decay becomes more efficient. On the other 
hand, interband decay becomes less probable owing to 
the reduced overlap with bulk states. It turns out that 
the increase of intraband decay outweighs the decrease 
in interband transitions. The net result is a more steep 
increase of T than obtained from parabolic dispersions. 
We thus find the intraband contribution to dominate up 
to E — Ep = 1 eV where the interband and intraband 
contributions become similar in magnitude. Comparing 
to published experimental data (symbols in Fig. [3Jl the 
agreement is improved. 

In summary, we have presented a theoretical many- 
body analysis of the inelastic e-e lifetimes of electronic 
excitations in the Shockley surface state band of Ag(lll). 
The calculations are based on the GW approximation us- 



ing a one-dimensional pseudopotential together with an 
ab-initio description of the projected surface band struc- 
ture. We have shown that the consistency of experimen- 
tal and theoretical derived damping rates is significantly 
improved when using ab-initio band structures instead of 
an effective mass approximation for calculating the decay 
rate. 
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